Exciton BCS or BEC state in a semiconductor bilayer system? 
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We calculate the off-diagonal long range order (ODLRO) terms of the exciton-exciton correlation 
function of a semiconductor bilayer system with Coulomb interaction and a transverse magnetic 
field. We show that the formation of a BEC state is very sensitive to the width of the interaction in 
momentum space. This dependence is analytically derived and represents the key physical ingredient 
for the formation (or not) of an exciton condensate state. 
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INTRODUCTION 



LANDAU LEVELS 



After the formulation of the pairing model for super- 
conductivity, speculations about the formation of exciton 
condensates startedfH- Much theoretical and experimen- 
tal effort was invested in the search of this condensate. 
However due to the short life-time of optically generated 
excitons the detection is far from being a simple task. 
A new approach to this problem has been then adopted 
with indirect excitons in semiconductor bilayer systems 
@) S 0] and it has been claimed Q that there is a BEC 
phase in the system. In this last experiment two parallel 
semiconductor plates with an applied external transverse 
magnetic field are under investigation. The presence of 
a magnetic field creates additional energy levels in the 
system, known as Landau levels, which have a degener- 
acy proportional to the magnetic flux piercing the plates. 
Therefore if the magnetic field is strong enough, there 
will be more available lowest Landau levels (LLL) than 
the number of electrons, such that the unoccupied levels 
can be seen as holes in resemblance to the usual holes in 
semiconductor valence bands. 



Many theoretical papers have dealt with this kind of 
system as in Refs. p, 0, 0] to cite a few. Nevertheless, a 
theoretical approach to such system is still a great chal- 
lenge, as it is for any many-body problem. It is too hard, 
if not impossible, to determine the ground state of such 
systems and therefore either approximations or numer- 
ical calculations are called for. In the case of an exact 
numerical treatment of the GS wave function this is re- 
stricted to a small number of particles. Here we take 
a look at the experiment described above and propose 
simple models for the ground state of that system which 
allow analytical results and point to the key physical in- 
gredient responsible for the BEC formation. We calculate 
the off-diagonal long range correlations for the excitons 
since it characterizes the formation of a condensate. 



Consider an electron moving on the x-y plane under 
a constant magnetic field B in z direction. One possible 
gauge for the potential vector A is 
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Here we find that V • A = 0, which implies that p and 
A commute. Writing down the Schrodinger equation, we 
find 
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ip(x,y) = Etp(x,y) 



We try a solution of the type 
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which yields 
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Xn(y) =E nXn (y) (4) 



where oj c = eB/(mc) is the cyclotron frequency and I = 
y/nc/JeB) is the magnetic length. The solutions are the 
same of a harmonic oscillator centered at yo = —l 2 k and 
are given by 



X n (y)=e-(y +l2k ?I^H n [{y + i 2 k)ll] 
V 
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Considering only the LLL, i.e., n = 0, the normalized 
wave-function is 



ipk(x,y) 



1 
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These wave-functions are orthogonal and obey 

d 2 r ijjl(r)i>k' (r) = Sk,k> , 
Y^r k {r)Mr')=5(r-r') . 
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In a finite system, k runs from 1 to Q, which is the degen- 
eracy of the system. In the present situation of a plate of 
area A — L x L y , we find £1 = BA/$ , where $o = hc/e 
is the elementary magnetic flux. 



Second quantization formalism 

Using second quantization, we define the field opera- 
tors for the electrons as 

*(r)^^.(r)a k , *t (r ) = £ ^ (r)a t ; (9) 



where a k (at) are the annihilation (creation) operators 
in the Fock space. The action of these operators in a 
Fock state is 

akK ■■■n k ---n n ) = n k \m ■ ■ - 0- • -n n ) , (10a) 
a k |ni • • -n k ■ ■ 'Tin) = (1 - n k )\ni -no) , (10b) 

and they obey the canonical anticommutation relations 

{a k ,a k ,} = S k k' , {a k ,a k /}=0, {a k ,a k ,}=0. (11) 

The total number of particles in such a state is obviously 
./V = y\ rtfc and this can be written in an operator fash- 
ion as 



N = ^a k a k = /d 2 r* f (r)*(i 



(12) 



Let us compute the Fourier transform of the density 
operator p(r) = ^(r)^(r) for electrons given by 



p(q) = J d 2 rc- lq r * t (r)*(r) 

= n l/2 L l E j d ' 
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In the last passage we have defined Ak = k 1 — k and 
k = (k + k')/2. Integrating we get: 

p(q) =J2^,Ak e- i2Afc2 / 4 exp{-^ + z 9? /fc}aW 
= e - ;2 9 2 /4 ^ qyq j2 e u2 k *y a t k a k+qx 

k 

(14) 



Particle— Hole transformation 

We are initially dealing with electrons in two parallel 
sheets. However if the number of electrons is smaller than 
f2, we can treat the empty levels in one of the sheets (the 
lower one for instance) as holes. We then define the hole 
creation operator ^(r) as the electron destruction one. 
Thus 



\&t(r) = *(r) and h\ = a k 



(15) 



Then the Fourier transform of the density operator for 
holes is similarly 



p(q) = J d 2 re^ q r * t (r)*(r) 



(16) 



COULOMB INTERACTION 

We are mainly interested in the Coulomb interaction 
between the particles in different sheets. The Fourier 
transform of the Coulomb potential in 2-D and a separa- 
tion in z equal d reads 



47re 2 e~ qd 
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The coulombic repulsion between electrons in the dif- 
ferent sheets can be more conveniently rewritten as a 
coulombic attraction between electrons and holes by a 
particle-hole transformation in one of the sheets. Then, 
this interaction term is given by 

- J d 2 r 1 d 2 r 2 * t (r 1 )*t( r2 ) 1 /( r2 _ n)*(r a )¥(rO 
= -(27r)- 2 ^p(q)U(q)p(-q) 



= Y, f p ^k^ic bic'+pak+p 
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where 



F = -- ?L- 

p ~ T c 



dq y cos(Z 2 q y q x ) e l2q2/2 



(19) 
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Motivated by the model proposed in Ref. [||, we first 
make a drastic approach to the interaction potential to 
be a delta function. From that, we find that F p oc 5 p q, 
and we have something similar to a pairing model. The 
ground state of this model is easy to find and we will use 
it to calculate the exciton-exciton correlation function. 
Later, we come back to discuss the real nature of the 
interaction potential and its implications. 



3 



GROUND STATE AND CORRELATION 
FUNCTION 

For the case that F p cx S p0 (pairing model) we can 
use the following expression to describe the ground state 
formed by n excitons: 



^(E4b£)"o>, 



(20) 



where a\, (b£) creates a electron (hole) in the up- 
per (lower) sheet and the normalization being J\f — 

m n \/(n-n)\. 

We now look for the general 2-body correlation func- 
tion as given by 

H H ^fei( r l) a L ^*k 2 ( r ^ h i 2 ^k' 2 (r' 2 )b k ' 2 ^ K {r' l )a. K \n) 



k2-k' 2 



(21) 



With use of Eq. (j2"0)) and the anticommutation rela- 
tions of the operators, we can calculate the expected 
value 



H a L b L h k' 2 a fc ; \n) = CiSk^h'^ + Ci^feife^*^ 



where 



C 1 = "ff~ n ; and C 2 = n(n - 1) 
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To see if this state is a BEC state, we calculate the 
exciton-exciton correlation function. This function is 
given by the creation of a pair electron-hole at one site 
and its following annihilation at another site, that is 



k,k' 



(24) 



fc,fe' 



First we note that V>fc(r) = ipl(r). Then there are 
only two products to calculate, namely, |V'fc( r )| 2 an d 
if>k{r)ipt(r'). We begin with the first one. Recalling 
Eq. 0, we find 



-l 2 (k + y/l 



2 ,2 
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We now transform this summation into an integral, with 
the measure dk = 2ir/L x . We find then 



^\^k(r 



1 
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For the second term we find 
1 



L x l^/n 



exp 
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Thus we get 

x|(0-n) + (n-l)exp(-^^)} (28) 

Taking the thermodynamic limit of the previous result, 
-> oo and n, — > oo, such that the filling 



with \r — 

factor v = n/fl remains finite, we get 



(n\W(r')W{r ')*(r)*(r)|n) 



(2ttZ 2 



(29) 



We then conclude that the ODLRO correlation terms 
remain finite in the thermodynamic limit, which charac- 
terizes a BEC state. That means, we would find a BEC 
if the interaction were sharp enough so that we could ap- 
proach the potential to a delta function. However, we 
know that the Coulomb potential is long range, and that 
changes the scenario completely. 



A MORE REALISTIC APPROACH 

Looking back to Eq. (fH)|) . by plotting F p for different 
values of d/l, we see that it has a typical width in p of 
the order of Ap ~ 1/ (10(H). This means that in this 
interval we find approximately VT2/250 values of p which 
contribute in the summation. In other words, this means 
that we will find excitons formed by an electron and a 
hole with different momenta, or simply that the exciton 
has a spatial 'width'. This is also ratified by a numerical 
solution for a system with a few states. 

With this in mind, we must generalize the state for the 
system in order to take this 'exciton width' into account. 
We then rewrite the state of n excitons with a separation 
in momentum p as 



1 



|0>. 



These states are still orthogonal, since 



(m,q\n,p) = 8 nm 5 pq 



(30) 



(31) 
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We can construct now a more general Ansatz for the 
state of the system as 



fi-i 

p=0 



(32) 



where ^2 p \f p \ 2 = 1 . We will assume for sake of simplicity 



that f p = e~P / 2q / v /(a v ^F+l)/2. 

To compute the exciton-exciton correlation function in 
this case, we first calculate the terms 



(p, n\^(r')^(r')^(r)^(r)\q, n) 



1 



(n — 1) exp 



2/ 2 
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+ (fi - n) jp&-*) e - ; v/2 U 
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Terms not diagonal in p, q yield a contribution of the 
order of I/O, which is negligible in the thermodynamic 
limit. The same happens to the first term of the above 
equation. We then have only to consider the following 
expression 

(*o|#V)*V)*( r )*( r )l*o) 
_ 1 n(n-n) ^ J P { X '- X ) 2e~P 2 / a2 Z 2 p 2 /2 



(2ttZ 2 ) 2 0(0 



a-y/TT + 1 



n(0 — n) 



(2ttZ 2 ) 2 0(0 - 1) V a 2 Z 2 + 2 



• exp 



2(a 2 l 2 + 2) 
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If a = 0, we reproduce the previous result of Eq. I2U1 
as we should expect, then a — is equivalent to make 



fp 



However, if a ^ 0, this term would vanish in 



the thermodynamic limit. 

From the previous result we must conclude that the 
condensate would be formed for the \n,p — 0) state. 
However, the solely Coulomb interaction is not sharp 
enough in momentum space to allow for a condensate 
formation. Based on the experiment related in Ref. [f|, 
we are led to conclude that either there is another in- 
teraction between the particles of the different sheets or 
the state they have found is not a BEC state, but just a 
collective BCS-like one. 
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